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Abstract 

The neutrino mass hierarchy, i.e., whether the neutrino mass eigenstate is heavier or lighter 
than the and z^2 mass eigenstates, is one of the remaining undetermined fundamental features 
of the neutrino Standard Model. Its determination would represent an important step in the 
formulation of the generalized model, and would have a profound impact on the quest of the 
nature of neutrinos (Dirac or Majorana) and the search for a theory of flavor. In this review, we 
summarize the status of experimental and theoretical work in this field and explore the future 
opportunities that emerge in light of the recently discovered non-zero and relatively large third 
neutrino mixing angle ^13. 
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1 Introduction 

Intense experimental and global analysis effort over the last several decades led to only a very small 
number of observations signifying that the standard model of elementary particles is incomplete and that 
ultimately it requires a generalization. Among them the existence of non-vanishing neutrino masses has 
a prominent position. Despite the obvious difficulty of experiments involving neutrinos, it is now clear 
beyond any doubt that neutrinos have masses, albeit tiny, and that they mix, unlike the assumptions 
of the standard model. Thanks to the discoveries of the last twenty years, we now have a simple 
three-flavor description of these new phenomena that explains most of the data. Despite the success, 
this is not the end of the story. Further experimental and theoretical effort is needed so that the next 
step, formulation of a new theory that encompasses the new insight into the neutrino physics, can be 
accomplished. 

Experiments with solar, atmospheric, reactor, and accelerator neutrinos have determined, with 
remarkable accuracy, most of the necessary parameters that describe the three-generation oscillations. 
However, several crucial pieces are still missing: the neutrino mass hierarchy (or the ordering of the 
neutrino masses), the magnitude of the CP (charge and parity) phase 6cp, the absolute scale of the 
neutrino mass, and whether neutrinos are Dirac or Majorana fermions. This review is devoted to the 
problem of the neutrino mass hierarchy. 

In the following, we assume that the neutrinos z/g, and Ur (so-called flavor eigenstates) that 
couple with the gauge bosons through weak interactions are coherent superpositions of three mass 
eigenstates z/j,i = 1,2,3, i.e. 

Ue2 U,3 \ / Z/i \ 

U,2 ■ Z /2 , ( 1 ) 

Ur2 U^3 J \^3 J 

where [/ is a unitary 3x3 matrix. In other words, we assume that there are three flavor and mass 
neutrino states and dismiss, for the purpose of this work, the possible existence of additional sterile 
neutrinos. This matrix implicitly contains possible charged lepton mixing; it is assumed that the 
corresponding charged lepton matrix matrix is diagonal in the flavor basis. 

The lepton mixing matrix t/ PQ |2l 13] (usually denoted as the PMNS matrix for Pontecorvo, Maki, 
Nakagawa and Sakata who were the pioneers of the field) can be parametrized in terms of the three 
angles ^ 13 , 623, and Ou and the CP phase 5cp- 

/ 1 0 0 \ / ci 3 0 \ / C 12 S 12 0 \ / 0 0 \ 

U = 0 C23 S23 ■ 0 1 0 ■ -S 12 C 12 0 ■ 0 0 

V 0 -S23 C23 / V 0 ci 3 /VO 0 1 / V 0 0 1 / 

C 13 C 12 C 13 S 12 Si 3 e"*'^crp \ / gia 0 0 \ 

-C23S12 - Si 3 S 23 Ci 2 e*'^C-P C23C12 - Si^S23Sl2p^^P C13S23 ■ j 0 6 *^ 0 ( 2 ) 

S 23 S 12 - Si3C23Ci2e*'^^-P -S 23 C 12 - Si^C23Si2p^^P C 13 C 23 ) \ ^ ^ ^ ) 

where e.g. Cij = cos 9ij and Sij = sin 9ij, etc. Here, a and (3 (both unknown at the present) are so-called 
Majorana phases that are decoupled from the phenomenon of neutrino oscillation. 
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Table 1: Neutrino oscillation parameters taken from Ref. |1]. For the atmospheric mass-squared splitting 
(Arug^,), the best fit results for both the normal and the inverted mass hierarchy are shown. These values 
are used in all the following plots except where noted. 


parameter 

best fit value ± la 

3a range 

sin^ 9 i 2 

u.oun_Q 0^2 

(0.270, 0.344) 

6*12 (degrees) 

33A8tou 

(31.30, 35.90) 

sin^ 023 

[0.4517^°}] or 0.5777^^^ 

(0.385, 0.644) 

6*23 (degrees) 

[42.2^;}] or 49.4^^;^ 

(38.4, 53.3) 

sin^ 6*13 

0.02197^°}? 

(0.0188, 0.0251) 

6*13 (degrees) 


(7.87, 9.11) 

6cp (degrees) 

251+®^ 

(0, 360) 

Am2i xl0“^ eV^ 


(7.03, 8.09) 

(normal) Aml^ xl0“^ eV^ 

+2.4587^^® 

(+2.325, +2.599) 

(inverted) Am|2 xl0“^ eV^ 

-2.4487®:®!^ 

(-2.590, -2.307) 


Table [T] summarizes the present knowledge of neutrino masses and mixings including neutrino mixing 
angles, the CP phase 6cp, and neutrino mass-squared differences Am]j = — m'j, based on the recent 

£t |lj after the Neutrino 2014 conference |5]. A comparable result has also been obtained in Ref. | 6 ]. 
The pattern of neutrino masses and mixings is schematically shown in Fig. At the present time, 
we cannot decide whether the neutrino mass eigenstate is heavier or lighter than the z/g and z /2 
neutrino mass eigenstates in nature. The scenario, in which the z /3 is heavier, is referred to as the 
normal mass hierarchy (NH). The other scenario, in which the z /3 is lighter, is referred to as the inverted 
mass hierarchy (IH). 

Superficially, it appears that the difference between the two mass orderings depicted in Fig. is 
striking, and hence it should be relatively easy to decide which one of them is realized in nature. As 
we show in detail further this is not the case and, in fact, the determination of the true mass hierarchy 
(MH) is very challenging. However, despite the difficulties, it is vitally important to solve this problem, 
and correspondingly a great effort is devoted to its solution. 

Lets hrst explain the nomenclature. Since the neutrinos has the largest component of the electron 
neutrino z/g while z /3 has the smallest component of z/g, the normal hierarchy in some crude sense 
resembles the mass ordering of the charged leptons, hence it is denoted as ‘normal’ . Obviously, the 
inverted hierarchy represents the opposite situation. Note that the word hierarchy is used throughout 
this work even if the lightest neutrino mass rriL (z/i for the NH and z/3 in the IH) could be comparable to 

the other two masses, i.e. when all three masses are much larger than the (so called degenerate 

scenario). 

Why is it important to determine the MH? The ultimate goal of all experiments and theories centered 
on neutrinos is to formulate the particle physics model that explains the observed neutrino masses and 
mixing patterns, and relates them to the well known charged lepton masses (and possibly to the quark 
masses and mixings). Clearly, explaining MH is a crucial part of that, and its determination would 
strengthen, or eliminate, roughly half of the proposed models. 

Neutrino masses can be also constrained or ultimately determined by cosmology and astrophysics, 
where the sum of the masses, rricosmo = ^mi is the relevant quantity. So far upper limits on that sum, 
more and more restrictive, were reported [7]. Depending on the MH, a lower limit of rricosmo is ~ 0.05 
eV for NH and ~ 0.1 eV for IH. Thus, knowing the MH changes the lower limit by a factor of two, and 
correspondingly restrict the range of the sum. 
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Most theoretical models of neutrino mass assume that neutrinos are massive Majorana fermions. The 
best way to test such a hypothesis is to search for the neutrinoless double beta decay Oi/jSjS] its rate is 
proportional to the square of the effective neutrino mass niee ■= \ \Uei\‘^mi + 

That quantity is restricted from below, ruee > 14 meV (taking into account the 3 ct error bars of the 
oscillation parameters) for IH while rUee = 0 is possible for NH. Thus, if IH is realized in nature, the 
next generation of the experiments can decide whether neutrino are Majorana fermions or not. 

Neutrino Mass Hierarchy 



Figure 1: Pattern of neutrino masses for the normal and inverted hierarchies is shown as mass squared. 
Flavor composition of the mass eigenstates as the function of the unknown CP phase 6cp is indicated. 
^^Itm ~ ~ l^rnhl and Am 2 i stands for the atmospheric and the solar mass-squared 

splitting, respectively. 


Similarly to most of the parameters describing neutrino mass and mixing, the neutrino MH can be 
accessed through the neutrino flavor oscillation. As shown in Table there are two small parameters in 
the neutrino oscillation description; the mixing angle ^13 (sin^ ^13 ~ 0 . 022 ) and the ratio Am^i/Amg^ 
(~ 3%). Due to this feature, most oscillation results are reasonably well described in the framework of 
mixing only two neutrinos, instead of three. In this case, the probability of flavor change in the vacuum 
and the oscillation length are given by 


F(z// ui>) 


sin^ 29 ■ sin^ 



Am^(eV^) ■ L(m) \ 
F(MeV) ) 


L 


vacuum 


(m) 


2.48 ■ F^(MeV) 
Am2(eV2) 


(3) 


and, obviously, the sign of Am^ (the mass hierarchy) cannot be determined in such case. 

Therefore, in order to determine MH, i.e. to hnd effects that are sensitive to the sign of Am|^ 
or Am| 2 , one has to either go beyond the vacuum oscillation or go beyond the simple framework 
of two-neutrino mixing. Correspondingly, there are two direct ways to determine MH. In the hrst 
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method, instead of the neutrino propagation in vacuum, one uses the propagation in matter where an 
additional phase from the interaction between neutrinos and the matter constituents makes the hierarchy 
determination possible. This effect is commonly referred to as the the Mikheyev-Smirnov-Wolfenstein 
(MSW) effect |Hl [HI [10] or the matter effect. In the second method, one explores the aforementioned 
small difference (Am^i) between and Am^ 2 - either case, carefully designed experiments are 

required. In the following, without going into details, we show the basic principles of both methods. 

The oscillation phenomenon has its origin in the phase difference between the coherent components of 
the neutrino flavor eigenstates caused by the finite masses rrii of the mass eigenstates Ui {i = 1,2,3). When 
neutrinos propagate in matter, an additional contribution to the phase appears. That phase affects 
the electron neutrinos and antineutrinos only, since only electron neutrinos/antineutrino interact with 
electrons in matter through the charged weak current (IT^-boson exchange). All three neutrino flavors 
interact with matter constituents through neutral current (Z-boson exchange) with same strength. 
Thus, these phases do not have effects on the neutrino oscillation. The corresponding effective potential 
and the matter oscillation length are 


W; = V2GfN„ 


L 


matter — 


27r 


1.7 X lO^(m) 
p{g/c-m^)Ye 


(4) 


where Gj? is the Fermi constant, N^. is the electron density in matter, p is the matter density, and 
Yf. ~ 0.5 is the electron fraction. 

For the illustrative case of two flavors of neutrinos propagating in matter, the equation of motion 
for neutrinos in the flavor basis then has the form 


. ^ f 

dx \ l^a 


2tt 


— K 


COS 26 I 1 

Lvacuum Lmatter 

sin 20 

^ 9T 

^i^vacuum 


sin 20 

^ 9 r 

^J-^vacuum 


0 


z/„ 


(5) 


where k, = sign(m 2 — m^). For antineutrinos the sign in front of l/Lmatter is reversed. Here one can 
see clearly that in matter, unlike in vacuum, the oscillation pattern depends on whether m 2 is larger 
or smaller than mi. Note that the sign of Am^i can be compensated by the sign of the cos20, i.e. by 
the octant of the mixing angle 9. The ordinary matter contains electrons, thus what really matters is 
whether is dominantly the lighter or the heavier one of the matter eigenstates (z/i or z/ 2 ). 

Therefore, from the study of solar neutrinos, where the matter effects are crucial, we know that 
Am^i > 0, provided that we define, as it is customary, that z/i is the mass eigenstate that is the 
dominant component of the electron neutrinos However, given the size of the matter oscillation 
length, it is clear that in order to use the matter oscillations in Earth to determine the sign of Am \2 
and very long baselines are needed. 

Matter effect in earth may be used for the determination of MH in two ways. The first method utilizes 
the appearance of new neutrino flavors with accelerator neutrino sources, in practice the appearance 
of the electron neutrinos in a beam that begins as the muon neutrinos The probability that a 
beam of propagating in the matter of a constant electron density W will be transformed into z/g as 
a function of the neutrino energy Ey and the neutrino travel distance L is 


Z/e) ~ Pq + PsinS + Pcos5 + -^ 3 ) 


( 6 ) 
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where 


Po 

Pz 

Ps'm 5 
Pcos (5 


n ^ sin^ 2013 . N « n 

^23p-3^sin P-1)A], 
cos^ 023 -—sin^(y4A), 


= —a 


A2 

sin 2012 sin 20 i 3 sin 2023 cos 0 i 3 sin Scp . 


= a- 


A{1 - A) 

sin 2012 sin 20 i 3 sin 2023 cos 0 i 3 cos Scp 


A{l-A) 


sin Asin(AA) sin[(l — A)A], 
cos Asin(AA) sin[(l — A)A], 


(7) 

( 8 ) 

(9) 

( 10 ) 


and where 


^ AmjiL 

4E. ’ 


^"^21 A AT 

Ot — —r - A — y^G pN^ 


Amil' 


When the neutrinos are replaced by antineutrinos the signs of the CP phase Scp and of the matter effect 
signifying the parameter A are reversed. Given known values of mixing angles, when measurements at 
optimally selected LjEy with both neutrinos and antineutrinos are performed, the matter effect and 
the CP violation can be disentangled and the sign of MH (i.e. the sign of A) can be determined. Fig. 
shows the oscillation probabilities vs. Ey at 1300 km baseline with distinguishing features of MH. A 
constant matter density is assumed. At 1300 km, the difference in oscillation probabilities between a 
constant matter density and realistic earth density profile is negligible. 

In addition, MH could be determined using the atmospheric neutrinos exposed to the matter effect 
by their propagation through Earth. Resonant oscillation HUDSIDSI occurs either for neutrinos in the 
case of NH, or antineutrinos for IH. A detector that could resolve v from v needs only to demonstrate for 
which state the resonance occurs. Detectors that can only determine neutrino flavor can take advantage 
of the difference in the atmospheric neutrino flux between v and v and the differences in interaction 
cross sections at the detection to determine the hierarchy. Details of the needed strategy and the 
corresponding challenges with these two approaches will be described in the next section. 

Beside the matter effect, MH can be determined by exploring the small difference (Am^^) between 
Am^^ and Am ‘^2 fh® three-flavor neutrino framework with neutrino and antineutrino disappearance. 
The survival probability for a neutrino of flavor I in vacuum is given by 


Pii — 1 — [A 21 sin^ A 21 -|- Ag 


31 


sin^ A 31 A 32 sin^ A 32 ] 


( 11 ) 


where A\j = 4 |f/^*p|f/b |2 


^ 1 ^ I 1 ^ I ......... —= AmjjL/iEy. This general formula, however, does not show 

explicitly the difference in the survival probability for the two hierarchies. In the case of practical 
importance, namely for the electron antineutrinos, the formula can be transformed into 


Fee = 1 — 2S33C33 — 4 c‘^ 13 Si 2 Ci 2 Au^ A21 -|- COs{ 2 A'i 2 ± 0ee) , (12) 

where the angle 0 that characterizes the difference of the two hierarchies is 


Sm 0ee = 


C 32 sin 2 A 


21 


y/l - 4sf2C?2 sin^ A 21 


, cos (pee 


-12 


COS 2 A 21 -|- S 


’12 


yi - 4 s52c52 W A 21 


(13) 


Here the ± sign in the last term in Eq. (12) contains the information regarding the mass hierarchy. The 


plus (minus) sign corresponds to the normal (inverted) mass hierarchy. The analogous formula for the 
survival probability of the muon neutrinos is 

F^^ = i-2|F^3r-(if/^ir+if/^2r)-4|F^inF^2rsin2A2i 

- 4|F^3p • (^^|F^i|^ + |F^2|" + 2|F^i|2|F^2Pcos2A2i^ • cos(2A32 ± 0^;.) , (14) 
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Figure 2: Electron neutrino and antineutrino appearance oscillation probabilities as a function of the 
neutrino energy ai L = 1300 km for the indicated four values of the CP phase 5cp- The matter 
density p and the electron fraction are assumed to be 2.7 g/cw? and 0.5, respectively, sin^ ^ 13 , sin^ 6 * 12 , 
sin^ 6 * 23 , and are assumed to be 0.0219, 0.304, 0.5, 2.4x10“^ eV^, and 7.65x10“® eV^, 

respectively. Top (bottom) two panels correspond to the normal (inverted) mass hierarchy. Left (right) 
panels correspond to the neutrino (antineutrino) oscillation. In the normal (inverted) mass hierarchy, 
the neutrino (antineutrino) appearance are enhanced, and the antineutrino (neutrino) appearance are 
suppressed. 

where U^i are the elements of the second row of the PMNS matrix and now the angle (j) is 

__ |^/.ipsin2A2i _ ^_ I^MiPcos2A2i + |?7^2p _ 

""" V\U,i\^ + \Uf.2\^ + COS2A21 ’ ^\U,i\^ + \U^2\^ + 2\U^,\^\U,2\^cos2A2i 

(15) 

By an accurate measurement of the survival probability and its L/E^ dependence, combining with 
the accurate knowledge of at least some of the parameters involved, MH (the sign in front of 0) can be 
determined. Detailed description of the issues involved will be given again in the next section. 

Furthermore, there is a possibility of determining MH based on the observation of core-collapse 
supernovae. There are at best only a handful of the galactic core-collapse supernovae per century and 
their neutrino emission and detection is a complicated matter. On the other hand, neutrino sky coverage 
by the sufficiently large detectors has been on for about three decades already, and will continue in the 
foreseeable future. There is, therefore, a good chance that a high statistics SN neutrino sample will 
eventually become a reality and that the neutrino mass hierarchy can be deduced from its analysis. It is 
possible, perhaps even likely, that the SN neutrino spectra are not universal; they might depend on the 
progenitor mass, accretion or rotation. Nevertheless, some general features exist. Only recently it has 
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been recognized that in the deep core region where the neutrino densities are very large the neutrino 
self-interaction could result in a hierarchy dependent collective flavor conversion. (For a review see 
Ref. [E]). Generally, this flavor conversion is not expected to affect the v and V spectra in the case of 
NH, but will alter signihcantly and differently the v and v spectra in the case of IH. In addition, the 
neutrino spectra will be affected by the resonant matter effect in the outer lower mass density layers of 
the supernova. Resonance occurs for neutrinos in the case of NH and antineutrinos in the case of IH. 
Also, the detected signal could be possibly altered by the matter effect during the transitions through 
earth (for general references see the review in Ref. [15]). Given the fast development of the supernova 
explosion and neutrino production and propagation theory, we will not describe the possibilities of MH 
determination based on the future SN neutrino observation in further detail. 

mcosmo=mi+m,+m 3 n\ =y | mf +1pm|+1 





Figure 3: The top left, top right, and bottom left panels show the rricosmos (fhe summation of neutrino 
masses), (the effective neutrino mass probed in the beta decay experiments), and rUee (the effective 
neutrino mass probed in experiments searching for the neutrinoless double beta decay) vs. the lightest 
neutrino mass rriL for both the normal and inverted mass hierarchy scenarios. The bands cover the 
allowed region given Sa neutrino mixing range in Table [T] and all possible values of Majorana phases. 


Finally, some limited information on MH can come from cosmology, the study of the beta decay 
endpoint, and the study of the neutrinoless decay. In cosmology, the sum of the neutrino masses 
rricosmos = '^irrii (top left panel of Fig. 


‘^cosmos 


rriL + -|- Am^i -|- 

rriL + \/m\ + |Am| 2 | + \/m\ + |Am|i 


(Normal Hierarchy) 
(Inverted Hierarchy) 


(16) 










































could be determined. If that sum is below the minimum mass of ~ 0.1 eV corresponding to the inverted 
hierarchy, the existence of the normal hierarchy is indicated. Similarly, in the study of the beta decay 
spectrum near the endpoint, if the effective neutrino mass rriy^ = ^J\Uel\‘^rn\ + |17e2|^ur2 + (top 

right panel of Fig. 


m 


Ve 


+ Ci^sl 2 ^'fn 2 i + (Normal Hierarchy) 

+ cf 3 cf 2 |Am|]^| + (Inverted Hierarchy) 


(17) 


is measured to be smaller than ~ 0.05 eV corresponding to the minimum value of the inverted hierarchy, 
the existence of the normal hierarchy is suggested. In analogy, if the neutrinoless decay is reliably 
observed and the deduced neutrino Majorana mass ruee = \ \Uei\‘^mi + |f/e 2 pm 2 e^“* + |f/e 3 pm 3 e^^*| 


m 


2 

ee 


(^12^13^1^ + c^3Si2\/nr| + Am^i cos 2a + cos 2 / 3 )^+ 

s?3 V^L + sin 2/3)2 

(0^2^13 \/nri + \ rn\ + |Am|]^| cos 2a + Si^mi cos 2 / 3 )^+ 

(ci3'Si2\/^i + I^^31 1 sin 2a + sin 2 / 3 )^ 


(Normal Hierarchy) 


(18) 


(Inverted Hierarchy) 


with a and /3 being the Majorana phases (bottom left panel of Fig. is again below its minimal 
value corresponding to the inverted hierarchy, the normal hierarchy is implied. Note, however, that 
the proportionality between the double beta decay rate and is only one, albeit the simplest, of the 
possible ways that the decay could occur. 

The rest of this review is organized as follows. In the next section, the various experimental methods 
with sensitivity to the neutrino mass hierarchy are described in more detail. In Sec. the proposed 
experiments are described and discussed. Finally, in the last section the entire subject is briefly sum¬ 
marized. 


2 Various Methods to Determine the Mass Hierarchy 

Two known ways to determine the neutrino mass hierarchy (MH) were introduced in the previous 
section. The hrst one depends on the charge-current forward coherent scattering between (or z/g) 
and the electrons in Earth (matter effect). The second one explores the small difference between 
and Am| 2 - Both methods are now of practical significance due to the recently discovered non-zero, and 
relatively large, third neutrino mixing angle 6 * 13 . In this section, we describe the practical approaches to 
resolve MH that use available neutrino sources, namely i) the accelerator neutrinos, ii) the atmospheric 
neutrinos, and iii) the reactor antineutrinos. We also review the statistical interpretation of the mass 
hierarchy sensitivity, which is crucial for understanding the physics reach of various experiments to be 
discussed in the next section. 

2.1 The “Large” Third Neutrino Mixing Angle ^13 

In this section we briefly review the experimental discovery of the non-zero value of 013 ~ 8.4°, which 
opens doors to MH determination. 

Historically, the hrst attempt to determine the value of ^13 was done by the CHOOZ uniin] and 
Palo Verde [IH] experiments in late 1990s and early 2000s. Both experiments were reactor antineutrino 
experiments searching for oscillations of Ug at baselines of ~1 km with a single-detector conhguration. 
No oscillations were observed in either experiment and an upper limit of sin^ 2^13 < 0.12 was set 
at 90% conhdence level (C.L.) by CHOOZ. Almost 10 years later, in 2011, several hints suggested 
a non-zero ^ 13 . The hrst one was based on a tension [19] between the KamLAND disappearance 
measurement (a reactor antineutrino experiment at an average baseline of ~180 km) and the solar 
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neutrino measurements (e.g. ratio of v^. charge-current to neutral-current interactions from SNO). 

Subsequently, accelerator neutrino experiments MINOS [20] and T2K |2I] reported on their searches 
of to Vf. appearance oscillation that is sensitive to 6 * 13 . In particular, T2K disfavored the ^13 = 0 
hypothesis at 2.5 ct I2I1- In early 2012, the Double CHOOZ reactor antineutrino experiment reported 
that the 6*13 = 0 hypothesis was disfavored at l.Ocr with one detector [ 22 ] , 

In March 2012, the Daya Bay reactor antineutrino experiment reported the discovery of non-zero ^13 
with a >5a signihcance [23] • About one month later, RENO conhrmed the Daya Bay discovery with a 
4.9cr signihcance [23] • And later in that year, Daya Bay increased the signihcance to 7.7a with a larger 
data set [25]. By then, non-zero ^13 was hrmly established. The left panel of Fig. shows the current 
global status of sin^ 2^13 measurements compiled with the latest results from each experiment. The 
precision of this important parameter is still being improved by the current-generation experiments. 
In particular, the right panel of Fig. shows the expected uncertainty of sin^ 20 i 3 from the Daya Bay 
experiment; by the end of the experiment it is expected to reach a better than 3% measurement precision 
of sin^ 2613 . 
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Figure 4: (Left) The current global status of the ^13 measurements. Results are taken from Kam- 
LAND+SOLAR [T9|, MINOS [26], T2K [21], Double CHOOZ [21], RENO [29], Daya Bay (uH) [30], 
and Daya Bay [21]. Among them, T2K and MINOS have assumed 6 cp = 0. (Right) Uncertainty on 
the Daya Bay measurement of sin^ 2613 iii the precision era (FY14-FY17, or 2015-2018) under different 
assumptions from Ref. [32] . 

2.2 Mass Hierarchy from Accelerator Neutrinos Appearance 

The accelerator neutrino experiments provided and continue to provide crucial inputs to the development 
of the 3-flavor PMNS model. Lederman, Schwartz, and Steinberger already in 1960’s discovered the 
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muon flavor neutrinos [33] at the Brookhaven National Laboratory (BNL) and showed that the lepton 
flavor appears to be conserved in the weak interaction. The tau flavor neutrinos the second to last 
elementary particle in the Standard Model experimentally conhrmed, was discovered by the DONUT 
collaboration at the Fermi National Accelerator Laboratory (Fermilab) [33] in 2001. The K2K long- 
baseline experiment [35] confirmed the observation of the neutrino oscillation at energies comparable 
to those produced in the atmosphere. The OPERA experiment provided the evidence for to 
appearance [33] • The MINOS long-baseline experiment [37] provided the most precise measurement of 
the atmospheric mass-squared splitting Aml 2 - 

Very recently, the T2K experiment [3^ discovered for the first time ever the neutrino appearance 
oscillation (z/g appearance from z/^), at more than 5a statistical signihcance, and confirmed the “large” 
non-zero ^13 discovered in the Daya Bay experiment [3E]- The T2K experiment 1271 together with the 
reactor neutrino experiments in Refs. [3^ 130] ITT] provide an initial interesting hint of the leptonic CP 
violation. The most precise measurement of the atmospheric mixing angle 623 also comes from the 
T2K experiment. Besides the possible determination of MH, the accelerator neutrino experiments are 
currently the only known reliable way to measure CP phase Scp, as well as determining the 623 octant: 
whether the z /3 contains more muon neutrino z/^ or tau neutrino Ur- Furthermore, the measurement of 
the third neutrino mixing angle ^13 from accelerator experiments together with the measurement of ^13 
from reactor antineutrino experiments provides a nice unitarity test of the PMNS matrix [T3] . 


Absorber 


Positive Focusing 


Decay Pipe 



Figure 5: A schematic view of the conventional neutrino beamline. The magnetic held is toroidal 
and assumed to focus positive charged particle. The muon neutrinos and antineutrinos are mainly 
generated through the two-body decay of pions and kaons. The electron neutrinos and antineutrinos, one 
of the main backgrounds for the electron neutrino/antineutrino appearance measurement, are mainly 
generated through the three-body decay of kaons and muons. The dashed line illustrates the direction 
of very short-lived 7r°. 

Accelerator neutrinos are generated by a high-energy proton beam striking a nuclear target to 
produce pions and kaons, which in turn decay into neutrinos. A schematic view of a conventional 
neutrino beamline is shown in Fig. [^ First, energetic secondary pions and kaons are produced when 
high-energy proton beams interact with the nuclear target. Second, some of the charged pions and 
kaons within certain momentum range are focused by magnets (also known as horns), so that they 
are approximately traveling in parallel with the incident proton beam direction. The polarity of the 
magnet can be selected by changing the current in order to focus either the positive or negative charged 
particles. The charged pions and kaons then travel through a long decay pipe to provide enough time 
for them to decay. Finally, a thick absorber is placed at the end of decay pipe to absorb the muons 
(decay products of pions and kaons) and other remaining charged particles. The ” on-axis” neutrino 
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Figure 6: Neutrino fluxes for different species are shown for the Fermilab Booster neutrino beam¬ 
line [S]. Left (right) panel represents the positive (negative) polarity which focuses positive (negative) 
charged particles in the magnetic horn. Besides the intrinsic and fluxes, there are also sizable 
wrong-sign neutrinos (z/^ in the positive polarity mode and in the negative polarity mode), which are 
produced by the wrong-sign charged pions and kaons. Reproduced with the permission from Ref. [H]. 


beam refers to the case where the neutrino beam direction to detector is in parallel with the incident 
proton beam direction. In comparison, the ” off-axis” neutrino beam refers to the case when there is a 
angle difference between these two directions. The “off-axis” beams pioneered by the E889 proposal at 
BNL |32] result in a narrower neutrino energy distribution. 

The primary neutrino species are z/^ (or V^) through positive (negative) charged pion and kaon 
2-body decay: 


TT -F (vr —>■ p -F z/^) B.R. ~ 100% (19) 

+ {K~ —)■ /i“ -F Pfj) B.R. ~ 63%, (20) 

where B.R. stands for the branching ratio. Values are taken from Ref. [7]. There are also small amounts 
of irreducible and Pe, mainly produced in the 3-body muon and kaon decays: 

— y -F r'e “F z/^ P —y e -F r'e T B.R. ~ 100% (21) 

—>■ 7r° -F e"*" -F r'e 7 i° + e~ + Pe B.R. ~ 5% (22) 

+ e~ + Pe —)■ 7r“ -F e"*" -F z/e B.R. ~ 41%. (23) 

These intrinsic z/g and Pe in the beam are one of the main backgrounds for the z/^ —>■ z/g and z?^ —?• Pe 


appearance measurements. As an example. Fig. shows the neutrino fluxes at the Fermilab Booster 
neutrino beamline, where the incident beam proton kinetic energy is about 8 GeV. 

The accelerator-produced neutrinos are generally detected through the charged-current interaction 
with nucleon/nuclei in the detector. In an oscillation experiment, there is typically a detector placed 
near the neutrino production region which is intended to measure the essentially unoscillated neutrino 
flux and a second detector placed far from the source to be sensitive to oscillations. The commonly used 
detector technologies for far detectors include water Cherenkov, liquid scintillator, magnetized iron plus 
scintillator, and liquid argon time proiection chamber. The involved target nucleon/nuclei are typically 
^H, 160, ^OAr, and ^^Fe. 
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Fig. El shows the neutrino-nucleon and antineutrino-nucleon interaction cross sections. There are 
three main components based on the assumption that the neutrino interacts with one of the bound 
nucleons (“N”) in the target. 

In the quasi-elastic scattering (QE) 

z/ + ‘W” + N' 9+ “iV” ^ + N', (24) 

the nucleon is converted into its isospin partner. In the resonance scattering (RES) 

z/+‘W” ^+ iV* P+“iV” ^+ AT*, (25) 

the nucleon is excited to a resonance state. And in the deep inelastic scattering (DIS) 

V + “g” ^ /- + X P + “g” ^ /+ + X, (26) 

the neutrino is assumed to interact with valence or sea quarks (“g”) in the nucleon. 

Besides these reaction mechanisms, neutrinos can also interact with more than one of the bound- 
state nucleons (i.e. meson-exchange currents and short-range correlations). In addition, neutrinos can 
also interact with the entire nucleus by the elastic scattering 

z/ + A ^ + A' P + A ^ + A', (27) 

where the A’ is the daughter nucleus or by the coherent meson production 

v + A ^ l~ + A + v + A^l^ + A + rrT, (28) 

where m represents a meson (such as tt or p). The cross sections of the elastic and the coherent scattering 
in general are much smaller than the QE/RES/DIS cross sections. 




Figure 7: Summary of neutrino-nucleon interaction cross sections is reproduced with the permission 
from from Ref. |l5]. Left (right) panel shows the ratio of neutrino (anti-neutrino) cross section with 
respect to the neutrino energy. Three main reaction mechanisms based on the assumption that the 
neutrino interacts with one of the bound nucleon (“N”) in the target: quasi-elastic (QE), resonance 
(RES), and deep-inelastic scattering (DIS) are shown in addition to the total interaction cross section. 

For neutrino oscillation measurements, the knowledge of the neutrino energy E^, is crucial, as the 
neutrino oscillation probability is expressed as a function of and the neutrino traveling distance L. 
Depending on the reaction mechanism, there are different strategies to reconstruct the neutrino energy 
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(see e.g. Ref. [1^ and the references therein). For example, the qnasi-elastic scattering is assnmed 
to be a two body to two body scattering. For nncleons at rest in the lab frame, the nentrino energy 
can be reconstrncted from the charged lepton kinematics (momentnm and direction) together with the 
incident nentrino momentnm direction. For the target nncleons bonnd in the nnclens, the nentrino 
energy reconstrnction at QE kinematics z/; + “n" —)■ / + p can be written in this crnde bnt simple and 
often used approximation as: 


_ 2(M„ -Es)E,- [(M„ - + W - l^r-] 

•' 2[M„-Eb-E,+p,cos{0,)] ' ^ ’ 

with Mn, Mp, and Mi are the mass of neutron, proton, and charged lepton, respectively. Ei, pi, and 
9i are the energy, momentum, and polar angle of the scattered charged lepton. Eb is the separation 
energy. The above method has a systematic bias, as the reaction mechanism should include the Fermi 
motion, involvement of more than one nucleon in the initial state (such as the short-range correlation), 
and hnal state interactions. The correct modeling of the reaction mechanism in the neutrino event 
generator is thus crucial. Another example is the deep inelastic scattering. The neutrino energy can be 
reconstructed by the calorimetry method, which sums the charged lepton energy with the energy of the 
hadronic system. In this method, the correct modeling of the hadronic energy response in the simulation 
and the correct estimation of the missing energy due to neutral particles (neutron and neutrino) becomes 
essential. 

Fig. I^shows the calculated Up to z/g appearance oscillation probability vs. the ratio LjEy. Due to the 
matter effect, the different MH leads to different behavior for the neutrino and antineutrino oscillations. 
In particular, at the hrst oscillation maximum L/Ej, ~ 510 km/GeV, the normal (inverted) MH would 
enhance the neutrino (antineutrino) oscillations, but suppress the antineutrino (neutrino) oscillations. 
Therefore, by precision measurements of the Vp to z/g and the Vp to z/g oscillations, MH can be determined. 
Comparing to Fig. which is evaluated at hxed baseline L = 1300 km. Fig. |^is made by scanning L 
at hxed neutrino energy Ey = 2 GeV. Due to this choice, the matter effect is signihcant in the second 
and higher oscillation maxima in Fig. On the other hand, it is not as important for the second and 
higher oscillation maxima in Fig. as the neutrino energy is much lower. 

For neutrino oscillation experiments, the neutrino havor at detection can be tagged by the charged 
lepton havor. For the appearance experiments, the background determination is obviously essential. 
Besides the intrinsic Ve and z>g beam backgrounds, there are two additional backgrounds. The hrst one 
is the neutral-current 7r° production 


z/ -1- ‘W” ^ z/ + 7r° + X. (30) 

Due to its very short lifetime (~8.5xl0“^^s in its at-rest frame), tt" would immediately decay after 
production. About 99% of tt° would decay into two gammas. If one of the gammas is missed in the 

reconstruction (if it has low energy or because it is too close to the other gamma), the 71° could be 

misidentihed as an event. The second one is the charged-current scattering, when the neutrino 
energy is high enough to produce a r lepton, 

Vr + “A^” ^r- +X VrE “X” ^ r+ + X. (31) 

About 18% of would decay through 

r“ —)■ e“ -1- he + r’*' —>• -|- z^e + hr, (32) 

in which the e^ could mimic the z^e/hg charged-current interaction. 
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Figure 8 : The to z/g appearance oscillation probabilities are shown vs. the ratio L/E^. The matter 
density p and the electron fraction are assumed to be 2.7 g/cw? and 0.5, respectively. sin^ 6 *i 3 , 
sin^ 6 *i 2 , sin^ 6 ' 23 , Arug^, and Am^i are assumed to be 0.0219, 0.304, 0.5, 2.4x10“^ eV^, and 7.65x10“^ 
eV^, respectively. Top (bottom) panels correspond to the normal (inverted) mass hierarchy. Left (right) 
panels correspond to the neutrino (antineutrino) oscillation. The plot is made by scanning distance L 
at a hxed neutrino energy Ey= 2 GeV. Different curves represent different values of the CP phase 5cp- 


2.3 Mass Hierarchy from Atmospheric neutrinos 

In 1998, the Super Kamiokande collaboration reported a zenith angle dependent deficit of atmospheric 
muon neutrinos, inconsistent with the expected atmospheric neutrino flux This seminal result 

pioneered the entire held of the neutrino havor oscillations. Now, with the next-generation gigantic 
detectors, atmospheric neutrinos provide a new opportunity to determine MH. 

Atmospheric neutrinos are generated by high-energy cosmic rays colliding with nuclei in the upper 
atmosphere. The collision produces a shower of hadrons dominated by pions with a small amount of 
kaons. Similarly to the accelerator neutrinos, atmospheric neutrinos are generated through the decay 
of pions, muons, and kaons. The positively charged pions decay to a positively charged muon and a 
muon neutrino. The positively charged muons further decay to a positron, an electron neutrino, and a 
muon antineutrino. The decay chain of negatively charged pions is analogous. Therefore, one naively 
expects a two-to-one ratio of muon neutrinos to electron neutrinos, when the momentum distribution 
of decaying pions and muons is neglected. Fig. shows the comparison of the atmospheric neutrino 
flux calculations for the Kamioka site for different neutrino flavors, where above 2 GeV neutrino energy 
the ratio of muon neutrino to the electron neutrino begins to increase. 
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Figure 9: Atmospheric neutrino flux prediction is reproduced with the permission from Ref. m- 
The left panel shows the atmospheric neutrino fluxes calculated for the Kamioka site averaged over 
all directions. The right panel shows the corresponding flux ratios. Four models, HKKMll m, 
HKKMS06 |1H], Bartol [IH], and Fluka [HD] are compared. 


In addition, the ratios of neutrinos to antineutrinos are larger than unity. This is because the 
positively charged pions, containing the u and d quarks, are expected to have a higher production yield 
than the the negatively charged pions containing the d and u quarks. This, in turn, is caused by the 
fact that the primary cosmic rays are typically the high-energy protons, containing two u and one d 
valence quarks. 

The zenith angle of the atmospheric neutrinos is directly correlated with the neutrino traveling 
distance for a detector near the earth surface as shown in the left panel of Fig. The neutrino 


detection relies on the charged-current neutrino-nuclei interaction as discussed in the previous section. 
In particular, as shown in Fig. the neutrino interaction cross section is much larger (about factor of 
2 at high energies) than the antineutrino interaction cross section. 

Similarly to the accelerator neutrino case, MH can be determined through the matter effect as 
elaborated in Ref. [CT IS21 ESI ESI ES] among others. It is necessary, however, to take into account that 


the earth matter density is not constant, as shown in the right panel of Fig. 10 


Figs. and p^^show the calculated oscillation probability for the atmospheric muon neutrinos and 
electron neutrinos, respectively. In NH scenario, we have the following features 


• Muon neutrino disappearance and muon antineutrino disappearance there is 

a resonance for the neutrinos at energy b GeV and zenith angle cos^ ~ —0.95 (neutrinos 

traveling from the other side of the earth) compared to the antineutrinos. 


• Muon neutrino to electron neutrino appearance —)■ Ve and muon antineutrino to electron 

antineutrino appearance ^ Ve- there is a much larger oscillation in the neutrinos than in the 
antineutrinos. 


• Electron neutrino disappearance Ve —t Ve and electron antineutrino disappearance Vf. ^ ^e'- there 
is a much larger oscillation in the neutrinos than in the antineutrinos. 
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• Electron neutrino to muon neutrino appearance —>■ and electron antineutrno to muon an- 

tineutrino appearance z/g —)■ there is a much larger oscillation in the neutrinos than in the 

antineutrinos. 


In IH scenario, the behavior of neutrino oscillations are exchanged with those of antineutrino oscillations. 



Figure 10: Left panel: a schematic view of the atmospheric neutrinos generation and detection. The 
zenith angle 6 of the neutrinos are indicated. Right panel: The matter density of the Earth is shown as 
a function of the radius. The data is based on the Preliminary Reference Earth Model (PREM) [56]. 


The atmospheric neutrino flux contains four neutrino species (z/^, z/g, Pe) and the detectors are 
usually only sensitive to the neutrino flavor. Therefore, one cannot disentangle the muon neutrino dis¬ 
appearance from the electron neutrino to muon neutrino appearance. Similarly, one cannot disentangle 
the electron neutrino disappearance from the muon neutrino to electron neutrino appearance. Never¬ 
theless, they can still be sensitive to the matter effect, as the modihcation of neutrino oscillations due to 
the matter effect is large. In addition, as the neutrino flux is higher than the antineutrino flux and the 
neutrino-nuclei cross section is larger than the antineutrino-nuclei cross section, one could determine 
MH with a detector that is not sensitive to the lepton (in practice muon) charge. The sensitivity to 
MH can be further enhanced with a magnetized detector which is sensitive to the muon charge. 
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Figure 11: Atmospheric muon neutrino oscillations calculated with the nuosc package IS7I are shown 
as a function of the zenith angle cos6' and neutrino energy E^. Top (bottom) panels show the muon 
neutrino disappearance (electron neutrino appearance). Left (right) panels show the neutrino (antineu¬ 
trino) oscillations. The normal mass hierarchy and zero CP phase = 0 is assumed. For the inverted 
mass hierarchy, the neutrino and antineutrino oscillation patterns are exchanged. 
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Figure 12: Atmospheric electron neutrino oscillations calculated with the nuosc package IS7I are 
shown as a function of the zenith angle cos 9 and neutrino energy Ey. Top (bottom) panels show the 
electron neutrino disappearance (muon neutrino appearance). Left (right) panels show the neutrino 
(antineutrino) oscillations. The normal mass hierarchy and zero CP phase = 0 is assumed. For the 
inverted mass hierarchy, the neutrino and antineutrino oscillation patterns are exchanged. 

2.4 Mass Hierarchy from Reactor Antineutrinos 

Reactor antineutrino experiments have a glorious history. In 1956, Cowan and Reines first experimen¬ 
tally discovered neutrinos at the Savannah River reactor power plant in the U.S. [5H]. In 2005, the 
KamLAND experiment in Japan observed the neutrino oscillations in the solar sector [59]. Recently, 
the Daya Bay experiment in China, the RENO experiment in Korea, and the Double-Chooz experiment 
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in France reported the discovery of non-zero ^ 13 , the third neutrino mixing angle in 2012 |23l EH | 22 ] . In 
the future, besides the determination of the mass hierarchy, the reactor neutrino experiments will also 
provide precision measurements of neutrino mixing parameters including the mixing angles ^13 and Q \2 
and the mass-squared splittings Am ‘^2 ^^^1 

Nuclear power reactors produce electricity by the sustained nuclear chain reaction, and are essentially 
pure electron antineutrino z/g sources. In every nuclear fission ~200 MeV energy is released together 
with on average six z/g emitted by the /d-decay of the fission fragments. For each 1 gigawatt (GW) of the 
reactor thermal power, about 2 x 10 ^° z/g are emitted isotropically every second, making nuclear reactors 
one of the most powerful man-made neutrino sources, ^sspj^ together with the and 

produced by breeding in the reactor, are the main sources of the reactor z/g. The corresponding z/g 
energy spectra are shown in Fig. 
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^^®Pu |62], and 

^^^Pu [62] are shown. The curve with the positive slope represents the cross section of the inverse beta 
decay (IBD) process. The convoluted IBD spectrum, measured in experiments, is shown as the dotted 
line. 


The primary method to detect reactor z/g is using the inverse beta decay (IBD) reaction Ue+P —t e’''-l-n 
in which electron antineutrino interacts with a free proton. The final state particles are a positron and 
a neutron. An IBD event is represented by a coincident signal consisting of i) a prompt signal induced 
by the positron ionization and annihilation inside the detector such as a liquid scintillator (LS) detector 
and ii) a delayed signal produced by the neutron capture on a proton or a nucleus (such as Gd) with a 
large neutron capture cross section. In particular, the neutron capture on Gd releases multiple gammas 
with a total energy ~8 MeV. With 0.1% Gd-doped LS, the average time between the prompt and the 
delayed signal is about 30 ps. In comparison, for neutron capture on hydrogen a single gamma with 
energy of ~2.2 MeV is emitted. For liquid scintillator CnHr~. 2 n, the average time between the prompt 
and the delayed signal is about 200 ps. 

Fig. 14 illustrates the reaction and detection principle of the IBD process. The cross section of 


the IBD process is shown as the solid black line in Fig. lA with units showing on the right axis. 
The convoluted energy spectrum, i.e. the product of the reactor antineutrino flux and the IBD cross 


^There is also a small component of the electron neutrinos Vg with energy ^0.1 MeV from the neutron activation of 
the shielding materials. 
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section, is plotted as the dashed line. It begins at the threshold of ~1.8 MeV neutrino energy 
and peaks around 4 MeV. The tail of the spectrum extends beyond 8 MeV. Due to the time and space 
correlations, IBD can be easily distinguished from radioactive backgrounds which mostly consist of only 
a single signal. Furthermore, the energy of the prompt signal (also known as the visible energy E^is) is 
directly related to the neutrino energy: 

Ey = Ee+ +mn + Tn-mp^ Eprompt + 0.78 McV, (33) 

where mg, m„, and ntp are the mass of electron, neutron, and proton, respectively. is the relatively 
small (~ a few keV) kinetic energy of the recoiled neutron. Besides the energy of the positron, Eprompt 
also contains the mass of electron that annihilates with the positron. The straight-forward method of 
reconstructing the neutrino energy is a very attractive feature for the study of neutrino oscillations that 
require good knowledge of the neutrino energy. 


Scintillator 



Figure 14: 
from Ref. 


The detecting principle of IBD is shown. This figure is reproduced with the permission 


The basis for the MH determination using the hg disappearance has been extensively discussed in 
Refs. [6ll EH ESI EH EH ESI [TOl [TH [721 |73l [71], and is illustrated in Eq. (12) in Sec. In order to 
further elaborate on this point, we dehne the effective mass-squared splitting as: 


^^eff - T 


(2A32 ± 0 eg) — 


2 Am 32 -|- Am? 


2Am‘^2 


Ami 


(Normal Hierarchy) 
(Inverted Hierarchy) 


(34) 


which corresponds to the fourth term in Eq. (|12|) 
In the left panel of Fig 


well as the baseline L 
,2 
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Here, we define Am? := 4 ■ 0eg ■ EyjL. 


the Ami is plotted as a function of the visible energy E^is = Eprompt as 
At a Daseline L ~ 60 km, Am^ has a clear energy dependence. In particular, the 
Am^ at low energy (2-4 MeV) is larger than that at high energy (4-8 MeV), providing an opportunity 
to determine MH. For the normal (inverted) MH, the Am^jj measured at low energy (2-4 MeV) would 
be higher (lower) than that measured at high energy (4-8 MeV). However, in order to perform such a 
measurement successfully, the energy resolution of the detector must be very good (better than 1.9% 
for 5E/E at 2.5 MeV), since the oscillations corresponding to Am^^j are very fast at low energy and a 
cruder energy resolution would make the oscillation pattern to disappear. On the other hand, at short 
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baseline (L < 20 km), Am^ is essentially a constant, thus it would be impossible to determine MH 
through comparing values at low and high energy regions. The purple line in the left panel of 

Fig. lA represents the approximate boundary of degenerate mass-squared difference. The right side of 
the purple line alone will yield negligible contributions to MH. 

Alternatively, if one could measure Am^jj very precisely in both the electron antineutrino disappear¬ 
ance channel and in the muon neutrino and antineutrino disappearance channels, one can be sensitive 
to MH [75l |66]. This is because the Am^ values in these two channels are different. More specihcally, 
as shown in the right panel of Fig. 
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Am^ ee larger than Am^ 


Therefore, we have 


^^Iff ee = 2Ami2 + Ami ee > ^^Iff + Ami MM 


ee = 2Am22 - Am^ < Amlfj = 2 Am‘l^ - Ami mm 


(Normal Hierarchy) 
(Inverted Hierarchy) 


(35) 


Since the difference between the Am^ in these two channels is quite small (~ 5 x 10“® eV^), the 
required precision of Am^jj in both channels would be about ~ 1 x 10“® eV^, which is very challenging. 
Therefore, this method is likely to serve as a validation once the MH is already determined by other 
means. 



Figure 15: (Left) The effective mass-squared difference Am^, as a function of baseline (y-axis) and 
visible energy ~ — 0.8 MeV (x-axis). The color represents the size of Am^ in eV^. The solid, 

dashed, and dotted lines are three choices of detector energy resolution with 2.8%, 5.0%, and 7.0% at 1 
MeV, respectively. The left side of these lines (low values of E^is) will yield negligible contributions to the 
differentiation of the MH due to large values of LjEy The purple solid line represents the approximate 
boundary of the degenerate mass-squared difference. The right side of the purple line alone will also yield 
negligible contributions to the differentiation of the MH. (Right) Am^ for both electron antineutrino 
{Ami gg) and muon disappearance {Ami ^^) are shown as a function of the unknown CP phase Scp- 
Reproduced with permission from Ref. m- 

2.5 Interpretation of Mass Hierarchy Sensitivity 

The determination of the neutrino MH is equivalent to hnding the sign of Am \2 = Thus, 

the determination of MH represents the test of two discrete hypotheses (normal hierarchy vs. inverted 
hierarchy or NH vs. IH in short). In this context, there are two different but related questions: i) given 
the data from an experiment, how to quantify the finding regarding MH? and ii) when designing an 
experiment, how to evaluate its sensitivity to MH? 
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1 - Probability 




Figure 16: (Left panel) The probability distribution function of AT with AT 11 is shown. The 
overlap is due to natural spread from statistical and systematic uncertainties. (Right panel) Several 
sensitivity metrics in the Bayesian statistics with respect to AT ^Xmh Eire shown. See the text for 
details. Reproduced with the permission from Ref. EH- 


The usual approach to the first question is as follows. Given the result of a measurement, we compare 
the data (represented by “x”) with expectations from both the NH and IH hypotheses. In practice, one 
commonly defines a test statistic ■= —2 log Lath, where is the likelihood of NH for data “x”. 
Tnh is further minimized over all nuisance parameters, including unknown parameters and systematic 
uncertainties, to obtain Similarly, another test statistic Tf^” is calculated from the likelihood of 

IH for data ”x”. In order to evaluate whether the experimental data favors NH or IH, a test statistics 
AT = T]^” — T(J*|p is defined. It is easy to see that a positive AT would favor the NH hypothesis, 
and a negative AT would favor the IH hypothesis. In addition, the absolute size of AT contains the 
information of how much data favor one hypothesis relative to the other one. 

Regarding the second question, when designing an experiment, the community sometimes report 
the average expectation AT which is calculated using the following procedure: 


• Assuming that the true MH is NH, a data set can be generated based on the best known exper¬ 
imental parameters. Such a data set does not include any statistical fluctuations and is usually 
called the typical data set or the Asimov data set eg. 

• When this typical data set is compared with the NH hypothesis, we naturally have T^|p = 0. 
Therefore, ATmh=nh = T^”. The subscript labels represent the assumption that the true MH 
is NH. Similarly, one can calculate ATmh=ih assuming that the true MH is IH. Here, AT can be 
interpreted as either the most probable value of AT or the average value of AT. 


In reality, the data from any experiments always contain statistical fluctuations as well as variations 
due to systematics. Therefore, if the experiment would be repeated many times, a distribution of AT 
appears. Due to the nature of two discrete hypotheses of testing MH, when the observed counts in the 
experiment are large enough, one can show that the distribution of AT approximately follows a Gaussian 


distribution with the mean and the standard deviation as AT, and 2'^|AT|, respectively [771 EHl 179] . 
For example, the left panel of Fig. IT shows the expected distribution of AT for experiments with 

11. There will be about 4.8% of experiments which favor a wrong MH 


\ATmr= 


MH=NH\ 


\ATmh= 


IH\ 


^The notation of or sometimes Ay^ is commonly used to represent AT. 
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Ar AT 

1 

4 

9 

16 

25 

P(IH;x) vs. P(NH;x) 

38% vs. 
62% 

12% vs. 
88% 

1.1% vs. 
98.9% 

0.034% vs. 
99.966% 

3.7 X 10"® vs. 
100% 

P(IH;x)/P(NH;x) 

0.61 

0.136 

0.011 

3.4x10“^ 

3.7 X 10-6 


Table 2: Probability ratios are shown with respect to Ar values. 


preference. Therefore even if AT is large, there remains a probability that a small AT is observed in 
an experiment. For example, in a situation where NH and IH are such that |AT| = 25 , the chance is 
about 1% that the experiment generates observed data favoring the wrong mass hierarchy. It should 
be noted that such a result naturally has a small AT value. 

2.5.1 Bayesian Approach to the Mass Hierarchy Determination 


As illustrated above given data from a measurement one can calculate Tj^h ■ —2 log^L^n) and Tjh '■ 
—21og(T/H). Tnh and Tjh are further marginalized over all nuisance parameters including unknown 
parameters and systematic uncertainties ot the experiment to obtain and iju , 

respectively. Assuming that the prior information of NH vs. IH is 50% vs. 50%, i.e. that both 
possibilities are equally expected, the probability ratio of the IH vs. NH, which is P(IH;x) vs. P(NH;x), 

1 1 1 X A —Ar /2 1 'xU A rrimarqinalized rrimarqinalized tHHi tt —IS.tI ‘2 • i 

can be calculated as e ' vs. 1, with Ar = ^ [ZZj. Here, e ‘ is commonly 

referred to as the Bayes factor. As illustrated in Ref. m, an approximation Ar AT can be made in 
practice. More generally, Ar can be explicitly calculated through MC simulations or other advanced 
integration techniques. Table lists a few values of Ar and their corresponding probability ratios. In 
the Bayesian approach, the final results are presented in terms of a probability ratio which is a natural 
and simple way to present results for the MH determination [771EO]. 

Ref. [77] also provides a few sensitivity metrics based on the Bayesian approach that illustrate the 
capability of an experiment; these are useful in designing experiments. One of them is Pmh=nh'- 


average probability to obtain the correct mass hierarchy. The second one is Pmh=nh^ which is based 
on probability intervals. The interpretation of Pmh=nh following: given that NH is the truth, 

90% of the potential data would yield a probability of NH, P(NH;data x), higher than Pmh=nh- Right 
panel of Fig. 16 shows the metrics Pmh=nh (dotted blue line) ar 
respect to AT, assuming AT \/S.Tmh=nh\ ~ \^Tmh=ih\ and Ar 


Pmh=nh (dotted red line) with 
AT. 


2.5.2 Prequentist Approach to the Mass Hierarchy Determination 


Ref. [7] (Particle Data Group) contains a very nice review regarding hypothesis tests. Given a null 
hypothesis Hq and the alternative hypothesis Hi, type-1 error rate a is defined as the probability of 
rejecting the null hypothesis Hq, if Hq is true. On the other hand, type-II error rate (3 is dehned as 
the probability of accepting the null hypothesis Hq, if Hi is true. Furthermore, according to Neyman- 
Pearson lemma, the likelihood ratio test is shown to be the most powerful test given that Hq and Hi 
are both simple hypotheses. This provides the basis to use AT dehned above to perform the MH 
hypothesis test. Before data are obtained, we choose the type-I error rate a, and this will decide the rule 
of the hypothesis test. The a value will determine the rejection region, and it will in turn determine the 
corresponding type-H error rate (3. In our MH problem, one can choose the null hypothesis to be NH, 

^In the marginalization process, one integrates the likelihood function over the entire parameter space of the nuisance 
parameters. 

■^A simple hypothesis should not contain any undetermined parameters. 
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and the alternative hypothesis would be IH. The result of this test will tell us whether the NH would 
be rejected or not, given the pre-dehned rule. Similarly, one should also perform a second hypothesis 
test by choosing null hypothesis to be IH with the alternative hypothesis being NH. The result of the 
second test will tell us whether the IH would be rejected or not. Unlike the Bayesian approach, the 
Frequentist approach does not directly address the question how much one MH hypothesis is favored 
over the other MH hypothesis given the experimental data. 

Furthermore, given experimental data and test statistics, one can obtain the p-value. The p-value 
represents the probability of obtaining potential data that are less compatible with the null hypothesis 
compared to the current experimental data. We emphasize that one should not confuse the p-value in 
Frequentist approach with the probability ratios in the Bayesian approach. For the MH determination, 
a nice feature exists for the p-value of the null hypothesis Hq being the disfavored MH hypothesis. In 
order to illustrate this point, we use the simple model shown in the left panel of Fig. Given AT from 


data and AT from the design of the experiment, the number of standard deviations which corresponds 
to the 1-sided p-value of the null hypothesis Hq being the disfavored MH hypothesis can be written as: 


{number of a) = 


(lArl + yjATI) 

2vAF 


> 


( 36 ) 


The equal sign is reached when |AT| = |AT|. In other words, given the observed data AT, there is 
a maximum p-value of the null hypothesis Hq being the disfavored MH hypothesis, which corresponds 
to in terms of the number of standard deviations. We should further point out that this 

maximum p-value does not depend on AT. 

The sensitivity of an experimental design can be represented using the type-I error rate a at certain 
type-H error rate (3. Given a design of an experiment, the type-I error rate a and type-H error rate (3 
are directly connected. In other word, a can be expressed as a function of /?. A smaller type-I error rate 

a would lead to a larger type-H error rate [3 and vice versa. The community commonly uses \J Ay;^ to 
represent the MH sensitivity of an experimental design. A recent discussion can be found in Ref. 


16 


the choice of \J Ay^ corresponds to choosing the 


In a simplihed model shown in the left panel of Fig. 
value of type-I error rate a as the sensitivity metric, while type-H error rate f3 is hxed at 0.5 {(3 = 0.5). 

Through investigating the statistical problem regarding Spin discrimination of the Higgs boson, 
whose nature is very similar to the MH discrimination, authors of Ref. [ H 2 ] suggested to consider the 
value of a{l3 = a), which is the value of a when a = (3, as the sensitivity metric. This can be compared 
to the conventional sensitivity metric a{f3 = 0.5), which is the value of a at /3 = 0.5. In the simplihed 


model as shown in the left panel of Fig. 16, the latter sensitivity coincides with the sensitivity metric 
Fmh=nh in Ref. [ 77 j , which is the probability of obtaining a AT favoring the correct MH. Right panel 
of Fig. 


16 


shows Fmh=nh with respect to AT (cyan dash-dotted line). 


In general, the high energy physics community is more familiar with the Frequentist approach 
than the Bayesian approach. In particular, there exist traditional criteria regarding the evidence and 
discovery of new physics: the simple 3a (evidence) and 5a (discovery) rule. If the 1-sided p-value of 
data when compared to the null hypothesis is smaller than 0.13%, one can claim that the evidence of 
new physics is observed. Similarly, if the 1-sided p-value is smaller than 3x10“^ , one can claim the 
discovery of new physics. For example, in the case of the non-zero 613 discovery [3H] , the null hypothesis 
is 6^13 = 0. In the case of the discovery of Higgs boson [ 8 ^ l 8 l] . the null hypothesis is the non-existence 
of Higgs boson. 

Recently, there have been some inspiring discussions [HSl ES] on whether the 5a rule should be 
universally applied as the criteria for discoveries. The 5a rule was invented to take into account i) look- 


^When transformed into number of standard deviations, a smaller p-value would give a bigger number. 
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elsewhere effect and ii) hidden unknown systematic uncertainties. For the MH determination, there is 
clearly no look-elsewhere effect. In this case, the discovery criteria should be adjusted accordingly. We 
further note that in practice it is strongly preferable that the designed experiment is not dominated by 
systematic uncertainties. 


3 Experiments 


In this section, we briefly describe the existing and planned experiments that have sensitivities to the 
neutrino mass hierarchy. In the hrst part, the current status is discussed. While the analysis of the 
atmospheric neutrino data by the Super-Kamiokande experiment slightly favors the normal hierarchy, 
the comparison of the experimentally determined Am |2 mass splitting using the muon and electron 
neutrino disappearance favors the inverted hierarchy. These two hints, that so far lack statistical 
signihcance, are obviously incompatible with each other. The running NOz/A experiment will be able 
to shed some light on this issue, but needs to collect more data. Also, there is a significant degeneracy 
between MH and the value of the CP phase 5cp due to the baseline of NOz/A. 

In the second part of the section we describe the proposed and planned experiments that promise 
to resolve MH. They are based on the study of the accelerator neutrinos (DUNE/LBNF), atmospheric 
neutrinos (PINGU, ORCA, INO, and Hyper-Kamiokande), and reactor neutrinos (JUNO, RENO-50). 
These next-generation experiments are in various stages of planning and funding. 



Figure 17: (Left) A schematic view of the Super-Kamiokande detector. Reproduced with the permission 
of Kamioka Observatory, ICRR (Institute for Cosmic Ray Research), the university of Tokyo. (Right) 
Preliminary result of Super-Kamiokande mass hierarchy measurement in terms of the negative-two-log- 
likelihood ratio (or — Xmin) plotted as a function of the unknown CP phase 6cp- The 

normal mass hierarchy is slightly favored by Ay;^ ~ 0.9, corresponding to about 60% vs. 40% probability 
with the Bayesian interpretation. This plot is reproduced with the permission from Ref. 0. 
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3.1 Current Status 


3.1.1 Atmospheric neutrino results from Super-Kamiokande 

The analysis of the atmospheric neutrino data from the Super-Kamiokande experiment currently pro¬ 
vides a hint of the mass hierarchy. The Super-Kamiokande detector (left of Fig. [ItI ) is a cylindrical 
50-kton (22.5-kton hducial volume) water Cherenkov detector located 1000 m underground in western 
Japan. The signal readout uses an array of regularly spaced large photo-multiplier tubes sensitive to 
single photo-electrons [HE]. From the pattern of the detected Cherenkov light, one can reconstruct the 
energy, the direction, as well as the particle type of the charged leptons produced by the charged-current 
neutrino-nuclear interaction. Using the principle illustrated in Sec. |2.3 Super-Kamiokande’s atmo¬ 
spheric neutrino data are sensitive to MH. Preliminary result was presented at Neutrino 2014 |83 IH9] 
and shown in the right panel of Fig. 17 The normal mass hierarchy is slightly favored by ~ 0.9. 
When interpreted in the Bayesian framework, this preliminary result corresponds to 60% (normal) vs. 
40% (inverted) probability in favor of the normal hierarchy. 


3.1.2 Global Analysis of the |Am| 2 | Mass-squared Splitting 


assuming Normal Mass Hierarchy 
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Figure 18: |Am| 2 l reported at the Neutrino 2014 conference are summarized assuming normal (left) 
and inverted (right) mass hierarchy scenarios. While Daya Bay measures the Vf. disappearance using the 
reactors, the rest of experiments measure the and/or disappearance using accelerator neutrinos 
(T2K and MINOS/MINOS-I-) and atmospheric neutrinos (Super-Kamiokande and IceCube). 


As illustrated in Sec. |2.4 MH information can be also extracted by comparing the measurement of 
the effective atmospheric mass-squared splitting, i.e. the value of Am 32 , between the electron and 
muon neutrino disappearance channels. For the normal (inverted) mass hierarchy, the effective mass- 
squared splitting in electron neutrino disappearance channel is larger (smaller) than that in the muon 
neutrino disappearance channel. In other word, the Am ‘^2 value extracted from these two channels 
would be more consistent with each other when testing the true MH hypothesis than the incorrect MH 
hypothesis. Figure 18 shows the extracted |Am| 2 | from Super Kamiokande [H7j in Japan, IceCube [10] 
at the South Pole, T2K [HI] in Japan, MINOS/MINOS-I- [92] in the U.S., and Daya Bay [EH ESj in 
China at the Neutrino 2014 conference under both normal and inverted MH scenarios. While Daya 
Bay measures the Ug disappearance from reactors, the rest of experiments measure the and/or 
disappearance. Figure lA shows that IH are slightly favored compared to NH. Assuming all results are 
uncorrelated, the fit of NH results yields a y^/AZJF value of 4.4/4. On the other hand, the fit of IH 
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results yield a x^/NDF value of 1.85/4. Therefore, the IH is favored by ~ 2.55, which corresponds 
to a probability ratio of 78% vs. 22% in the Bayesian framework. 


3.1.3 NOz/A 


The NOz/A (NuMI Off-axis Ug Appearance) experiment [M] is an accelerator neutrino experiment in the 
United States. It is currently taking data. NOz/A measures the z/g (z/g) appearance as well as the z/^ (P^) 
disappearance in a z/^ (P^) beam. The neutrino beam is generated by the 120 GeV proton beam from the 
Fermilab’s Main Injector. The system consists of two detectors, with the detector near the production 
point characterizing the neutrino beam and the far detector at the baseline of 810 km studying neutrino 
oscillations. Left panel of Fig. IT shows the overall configuration of the NOz/A experiment. The far 
detector is located at 14 mrad off-axis of the incident proton beam to obtain a narrow band neutrino 
beam around 2 GeV energy. The 810 km baseline corresponds to the maximal oscillation at 2 GeV. In 
addition, the long baseline offers the experiment a chance to determine the mass hierarchy through the 
matter effect as discussed in Sec. 12.21 



Figure 19: (Left) A map illustrating the NOz/A experiment. (Right) The detector conhguration of 
NOz/A from Ref. |95]. Reproduced with the permission from Ref. [96]. 


Right panel of Fig. 19 shows the detector configuration of NOz/A The near (0.3 kt) and far (14 
kt) detectors are structurally and functionally identical, despite the difference in their absolute sizes. 
Such a conhguration allows to minimize the systematic uncertainties in the neutrino beam, the neutrino 
interaction cross section, and the detecting efficiency. The near detector consists of 18 thousand PVG 
cells hlled with liquid scintillator and arranged in alternating X/Y layers with an overall dimension of 
14.3 X 4.2 X 4.2 m^. The light is transported through a wavelength shifting hber and then read out with 
avalanche photodiode. The far detector is approximately 60 times larger with the overall dimension of 
60 X 15.6 X 15.6 m^. Since it is located on the surface, the far detector is exposed to a large hux of 
cosmic muons. Nevertheless, a factor of 4x10^ suppression has been achieved to reject the cosmic ray 
background, taking into account the timing and the direction of the z/^ and z/g beam. Figure shows 
the expected MH sensitivity from NOz/A after collecting a total of 36x10^° protons on target (POT), 























which corresponds to about six years (three years neutrino and three years antineutrino running) of 
data taking at 750 kW beam power. The mass hierarchy resolution has a large dependence on the 
value of the unknown CP phase Scp- This is so because both the matter effect and the CP phase 
6cp can induce differences between neutrino and antineutrino appearance oscillations, which leads to a 
degeneracy between the MH determination and the determination of CP phase Scp for certain values 
of 6cp- At the most favorable CP values, the sensitivity of MH determination can reach about 3a. At 
this sensitivity, there is a 50% chance that NOuA will be able to reject the unfavored MH hypothesis 
by 3a. Meanwhile, 6.7% of measurements will yield a result favoring the wrong MH hypothesis. 


NOvA hierarchy resolution 



Figure 20: Mass hierarchy sensitivity from NOz/A (taken from Ref. |96j) is shown. The “significance 
of hierarchy resolution” is defined as y/ Ax^ ■= Vt as introduced in Sec. 


2.5 


3.2 Proposed and Planned Experiments 

3.2.1 DUNE/LBNF 

In this subsection the long baseline experiments that aim at the determination of the neutrino mass 
hierarchy, among other things, are described. There were two proposals, the Long-Baseline Neutrino 
Experiment (LBNE) [UT] in the US and the Large Apparatus for Grand Unihcation and Neutrino As¬ 
trophysics and Long Baseline Neutrino Oscillations (Laguna-LBNO or LBNO) [HH] in Europe, designed 
to separate the signal of the two hierarchies with matter effect. While the fundamental goal of the 
two proposals is similar, they differ in important details, in particular in the length of the baseline. 
Due to the limited available resource and expected costs of these projects, LBNE and LBNO collabo¬ 
rations combined efforts and merged into a new experiment located in the US: the Deep Underground 
Neutrino Experiment at the Long-Baseline Neutrino Facility (DUNE/LBNF). We concentrate on the 
DUNE experiment and will describe it in details. There is no firm schedule for DUNE at the present 
time. Tentatively, the construction of DUNE/LBNF might begin in 2017 and the neutrino beam might 
be ready by ~2023. 

DUNE is an accelerator neutrino experiment under development in the United States. It consists 
of a horn-produced broad band beam with 60-120 GeV protons at ~1 MW from Fermilab, >35 kton 
fiducial volume Liquid Argon Time Projection Chambers (LArTPC) 4850-feet underground at Sanford 
Laboratory in Lead, South Dakota, as the far detector, and high-resolution near detector(s) at Fermilab. 
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Figure shows the overall configuration of DUNE. In order to optimize the physics sensitivity to the 
CP phase 6cp, the baseline of DUNE is chosen to be 1300 km [21], which is longer than any existing 
accelerator neutrino experiments. The large underground LAr detector also represents an excellent 
opportunity for the underground science including proton decay searches, detection of atmospheric 
neutrinos and supernova neutrino burst. In particular, the detection of atmospheric neutrinos will 
provide additional information regarding MH. 



Figure 21; Illustration of the DUNE/LBNF (formerly known as the LBNE) experiment |1UU] . 


Normal Hierarchy 



Inverted Hierarchy 



Figure 22: The —)■ oscillation probabilities are shown for the normal (left) and the inverted (right) 
mass hierarchies. Different curves represent different values of currently unknown CP phase 6c p- The 
first and second oscillation maxima locate at about 2.5 and 0.8 GeV, respectively. The matter density 
p and the electron fraction are assumed to be 2.7 g/cm^ and 0.5, respectively. 

DUNE is designed to precisely measure parameters that govern —)■ z/g, z/^ —)■ z/g, —)■ z/^, and 

oscillations, including the second mixing angle 6 * 23 , the third mixing angle 6 * 13 , the CP phase 
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6cp, and MH. Figure shows the —>■ Ue oscillation probabilities for different values of 6cp for both 
normal and inverted mass hierarchies. The principle of the MH determination in DUNE relies on the 

! For NH, the z/g (Pg) appearance probability at first oscillation 
GeV) is enhanced (suppressed), while for IH, the opposite is true, i.e., the Ue ipe) 


matter effect, introduced in Sec. 
maximum (~2.5 

appearance probability at first oscillation maximum is suppressed (enhanced). At 1300 km, the matter 
effect in DUNE is large enough so that there is no fundamental degeneracy between the matter effect 
and the ^cp-induced CP effect. In other word, given enough statistics, these two effects can be totally 
disentangled. Furthermore, the precision measurements of the oscillation probability will also enable 
a measurement of the currently unknown CP phase 5cp that may hold the key to explain the large 
matter over anti-matter asymmetry in our universe cnn. 
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Figure 23: Event display of two neutrino events from Monte-Carlo. (Left) Charged-Current Quasi 
Elastic (CCQE) interaction with hnal states containing a muon and a proton. (Right) Charged-Current 
Deep-Inelastic-Scattering (CCDIS) with hnal states containing an electron shower, two protons (one 
being generated by neutron hnal state interaction), one positive pion, and one negative pion. 

The detector technology in DUNE is chosen to be the liquid argon time projection chamber (LArTPC) jlU2j . 
This technology was hrstly proposed by C. Rubbia for the neutrino physics [103] in 1977. As the most 
abundant noble gas in the atmosphere (~1.3% by weight), argon can be extracted economically, and is 
an ideal material for the TPC. Electrons have very high mobility in LAr |104] . At 500 V/cm electric 
held and 87 K temperature, the electron drift velocity is about 1.6 km/s |105111061110711108] . which 
is equivalent to about 0.64 ms for electrons to travel 1 m. Electrons have low dihusions in LAr |109j : 
the longitudinal and transverse dihusion coefficients are about 5 and 13 cm^/s |1101 llll[ I112[ 1113] . 
respectively at 500 V/cm electric held. For a 1.6 ms drift time, these numbers correspond to a position 
resolution (one standard deviation) of about 1.3 and 2.0 mm for longitudinal and transverse dihusions, 
respectively. LAr also has relatively high density (~1.4 g/mL at 87K), advantageous for neutrino 
experiments. 

LAr is also a good scintillation material. The scintillation light (emission peaks at ~128 nm mn) 
can provide the absolute event time to determine the drift distance along the electric held direction. 

The coordinate information perpendicular to the electron drift direction can be determined with a 
high resolution detector (e.g. wire planes). The spatial resolution could reach a few mm. Finally, the 
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ionization energy required to produce an electron-ion pair is about 23.6 eV [TM[TT3] . Some of the 
pairs will recombine in LAr, and the remaining electrons will be collected at wire planes. The charge 
collected per wire will be closely related to the energy deposition per unit distance (dE/dx), which can 
be used for the particle identihcation (PID). In addition, the event topology can be used for PID as well. 
LArTPCs are excellent detectors of neutrino interactions. Fig. shows event displays for two neutrino 
events from Monte Carlo as an example. The high-resolution nature of LArTPC is self-evident. 




Figure 24: The test power (p = 1 — /3 as introduced in Sec. 2.5), which represents the probability of 
accepting the correct (NH) hypothesis while excluding the incorrect (IH) hypothesis at the 3 and 5a 
level, shown as a function of exposure in kt-MW-years mg. The width of the bands represent the 
range of CP values. Sensitivities are for true NH (left) and true IH (right). 



Figure 25: The sensitivity for a typical experiment (blue line) is compared to the bands within which 
68 % (green) and 95% (yellow) of experiments are expected to fall |116j . The dashed lines represent 
the value of the \/AT metric for which the probability of determining the correct MH is 50% (cyan), 
98.9% (blue), or 3.7 x 10“® (black) in the Bayesian framework. The normal and inverted hierarchy are 
assumed to be true in the left and right plots, respectively. 


In order to suppress the uncertainties in neutrino interaction cross sections and the neutrino flux, 
a hne-grained neutrino detector will be installed at FNAL near the beam source. The near detector 
will measure the absolute neutrino flux via i) neutral-current -f e“ —?• -|- e~ elastic scattering, ii) 

inverse-muon decay —)■ /i“ + i>e for the high energy neutrino component, and iii) charged-current 
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quasi-elastic scattering + n ^ + p and + p ^ + n. In addition, it will also measure the 

spectra of all four neutrino species, z/^, Ue, and z/g through the Iow-z/q method |117l I118j to enable an 
accurate prediction for each species of the far/near flux ratio as a function of energy. The hne grained 
detector will be able to precisely measure the four-momenta of secondary hadrons, such as the charged 
and neutral mesons produced, and to constrain neutrino interaction cross sections, which are crucial 
inputs to the oscillation analysis. Furthermore, the precise measurement of these hadrons produced in 
the neutral- and charged-current interactions will signihcantly improve the background estimation in 
the oscillation signal. 

Fig. 24 shows the expected sensitivity of the MH determination as a function of the accumulated 
exposure in units of kt-MW-years. As introduced in Sec. |2.5.2[ the test power 1 — /? is plotted at values 


of type-I error a corresponding to 3 and 5a signihcance. Fig. 25 show the sensitivity in terms of test 
statistic AT as introduced in Sec. 2.5.1 At 300 kt-MW-years exposure, most of experiments will yield 
a MH determination better than 3.7 x 10“® vs. 1 in the Bayesian framework. 

LBNO, before merging into DUNE/LBNF with LBNE, was a proposed experiment in Europe. The 
initial neutrino beam is expected to be produced from CERN with a 400 GeV proton beam with a 
power of 750 kW. The primary choice of the far detector location is the Pyhasalmi mine in Finland. 
Various other sites are also under consideration. A multi-tens of kilotons double-phase liquid argon time 
projection chamber is expected to reside deep underground. The baseline from CERN to Pyhasalmi is 
2300 km, much longer than that of LBNE (1300 km). The z/^ —)■ z/g oscillation patterns are shown in 
Fig. 26 Comparing to those of 1300 km (Fig. 22), the matter effect (leading to the MH determination) 
at 2300 km is better due to longer baseline and corresponding higher neutrino energies. Therefore, the 
MH sensitivity is expected to be much larger than that of 1300 km. Given a similar overall exposure 
in kt-MW-years, it has a factor of ~ 4 increase in the test statistics AT. 
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Figure 26: The expected z/^ —)■ Ug oscillation patterns in LBNO at a baseline of 2300 km. The matter 
density p and the electron fraction Yg are assumed to be 3.2 g/cm^ and 0.5, respectively. 


3.2.2 PINGU, ORCA, INO and Hyper Kamiokande 


As explained in Sec. 2A, the study of atmospheric neutrinos with very large detectors may be used to 
determine the neutrino mass hierarchy. The long baseline of the upward going atmospheric neutrinos and 
the corresponding large matter effect, make this program possible. We briefly describe three proposed 
experiments that are aiming, among other things, to resolve MH. 
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The Precision IceCube Next Generation Upgrade (PINGU) is a proposed new multi megaton array 
for IceGube upgrade at the South Pole Station. Its primary purpose is the determination of neutrino 
mass hierarchy using the detection of atmospheric neutrinos. The initial design will have a dense array 
of 40 strings (20 m of average string spacing) with 60 optical modules each (5 m module spacing on 
the string), located in the Deep Gore region of the IceGube detector. This will allow reconstruction 
of the Gherenkov tracks of events, i.e. of the energy and the zenith angle, for a relatively low energy 
region, 5-15 GeV, where the effect of MH is most evident. Note that the zenith angle determination is 
equivalent to the neutrino path length determination. The initial design of the array, and description of 
the corresponding simulation can be found in Ref. |119] . The PINGU detector will also improve other 
measurements of neutrino oscillation parameters, in particular of the mixing angle 623 - It will be the 
world’s largest detector of the possible supernova neutrino burst. It will also be able to search for the 
indirect signal from the WIMP dark matter annihilation or decay in the Sun or Galactic Genter. In a 
favorable scenario, PINGU installation might be completed in 2021 or 2022. 

Although the atmospheric flux contains both neutrinos and antineutrinos that are difficult to sep¬ 
arate experimentally, as well as both and z/g, the large size of the PINGU detector, and the corre¬ 
spondingly large event statistics, makes it possible to overcome these challenges. At the relatively low 
energies, most of the secondary particles will be contained in the detector volume. Simulation indicates 
that ~ 90% of the charged-current z/g and induced processes, and of the neutral current events, can 
be identihed. 
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Figure 27: MH sensitivity of PINGU vs. the running time is shown. The blue line (tracks) is based 
on the identihcation of the muon track events and the green line is based on the hadronic and electron 
events. The red line represents the total sensitivity. For comparison the dashed red line show the ideal 
scaling without systematics. Reproduced with permission from Ref. mg. 
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As shown in Sec. 2.3, it tnrns out that the survival probabilities for a given and cosO^enith for 


the in the case of NH is essentially the same as the probability for with IH, and vice versa. 
This appears to be an insurmountable problem. However, the differences in the cross section and 
kinematics of the v and V interaction with nuclei, as well as the differences in the atmospheric v and 
V fluxes, makes it possible to separate the two hierarchies. As described in Ref. | 120 j . it is possible to 
identify the “distinguishability” region in the space of and cos 6 zenith where the effect of MH is most 
pronounced. It is these regions where the most useful information of the MH can be extracted. 


Simulation indicates that the PINGU sensitivity to the MH will reach a 3a signihcance with about 
3.5 years of data collection. In Fig. 27, the preliminary plot of the MH signihcance versus the PINGU 
detector running time is shown. Note that only when all available information is combined a signihcant 
determination of MH is possible. 


ORGA (Oscillation Research with Gosmics in the Abyss) |121] is a proposed option aiming to 
determine the MH using the deep-sea neutrino telescope technology developed for the KM3NeT project 
(A multi-km^ sized Neutrino Telescope) |122j . The experimental principle of ORGA is similar to that of 
PINGU. Instead of deep ice in the south pole, ORGA is expected to deploy large 3-dimensional arrays 
of photo-sensors to detect Gherenkov lights in the deep Mediterranean Sea. A 3-5(T MH sensitivity is 
targeted for a ~20 Mton-year overall exposure given the neutrino detecting threshold around 5 GeV. 


The India-based Neutrino Observatory (INO) |123l 1124] is a planned underground experiment at 
Pottipuram, Theni (India) to study the atmospheric neutrinos. The project has been approved by 
the Indian government with the tentative start of data taking by 2020. The proposed detector is 
a magnetized iron calorimeter with 52 kiloton mass and 1.5 Tesla magnetic held, divided into three 
detector modules. Each module contains 151 layers of 5.6 cm iron plates with the resistive plate 
chambers (RPG) as active detector elements in the gaps between the plates. INO will have excellent 
muon energy and direction resolutions as well as a charge reconstruction capability. Besides MH, INO 
will be also sensitive to 623 and Am‘^ 2 - Fig. 28 shows the projected MH sensitivity of INO vs. data 


taking time |125j . An improved analysis taking into account the hadron energy reconstruction together 
with all oscillation parameters marginalized, show an enhanced sensitivity m- The AT is expected 
to reach 9.5 with 10 years of data taking for sin^ 2^13 = 0.1 and sin^ 6'23 = 0.5 m- Note that the 
MH sensitivity strongly depends on the true values of 623 and 6*13 which, however, might be precisely 
known by the time the detector begins data taking. The MH sensitivity is signihcantly larger if the 
mixing angle 6*23 is in the second octant, i.e. if sin^ 623 > 0.5. The same also applies to MH sensitivity 
in PINGU. 
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Figure 28: MH sensitivities of INO with respect to the running time are shown as Axino •= Left 
and right panels assume the true MH is normal and inverted, respectively. Different lines represent 
different values of sin^ 6*23 and sin^ 26 * 13 . During the analysis, all the oscillation parameters are hxed at 
their input true values. Reproduced with permission from Ref. ra- 


Hyper-Kamiokande ina is a proposed experiment near the site of the previous Kamiokande and 
the current Super-Kamiokande detectors in Japan. The master plan at present suggests that the 
construction and the operation might begin in 2018 and 2025, respectively. The detector technology 
will be water Cherenkov similar as in Super-Kamiokande, only many times larger, with a total mass of 
1 Megatons (560 kiloton hducial mass). While the main experimental program at Hyper-Kamiokande is 
the precision measurement of —)■ and —)■ z/g oscillations with accelerator neutrinos at a relative 
short baseline (295 km) to search for CP violation in the lepton sector. Hyper Kamiokande is also 
sensitive to MH with atmospheric neutrinos. Fig. ^ shows the atmospheric neutrino MH sensitivity 
with ten years of data taking. Again, the MH sensitivity is signihcantly larger if the mixing angle ^23 
happens to be in the second octant. 
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Figure 29: The MH sensitivity of Hyper-Kamiokande with respect to the true value of sin^ 623 is shown 
as := AT. The blue band represents the uncertainty of the currently unknown CP phase 6 cp- A 
ten-years exposure is assumed. Reproduced with permission from Ref. m- 


3.2.3 Medium baseline reactor experiments: JUNO and RENO-50 


As described in the Sec. 2.4, reactor experiments at baselines ~ 60 km are well suited for the MH deter¬ 


mination. Two such experiments have been proposed and are currently at various stages of preparation 
and construction. Here we describe them in more details, using the JUNO experiment as the main 
example. 

Jiangmen Underground Neutrino Observatory (JUNO) is a reactor neutrino experiment under con¬ 
struction in Jiangmen City, Guangdong Province, China |128[ I129j . The groundbreaking for JUNO 
started in January 2015. The plan calls for the construction to last until 2017, assembly and instal¬ 
lation in 2018-2019 and start of the data taking in 2020. Fig. 30 shows the JUNO site. The JUNO 


detector will mostly receive Pg from two reactor complexes at Taishan and Yangjiang. The Taishan re¬ 
actor complex contains six reactor cores with a total thermal power of 17.4 GW. The Yangjiang reactor 
complex has two reactor cores with a total thermal power of 9.2 GW. There are two additional reactor 
cores (9.2 GW) planned at the Yangjiang site. The average baseline of JUNO is ~52.5 km with a RMS 
(root mean square) of 0.25 km. 

The central detector of JUNO is a 20 kton liquid scintillator (LS) detector with a total overburden of 
1850 meter water equivalent. Figure shows a conceptual design of the spherical LS central detector. 
The photo-cathode coverage is expected to reach 75%. Together with high performance LS (high photon 
yield with >14,000 photons per MeV, optical attenuation length of 30 m or better) and high quantum 
efficiency PMTs (peak quantum efficiency reaching 35%), JUNO aims to achieve an energy resolution 
of better than 1.9% at 2.5 MeV, which is essential for the MH determination. This level of energy 
resolution requires a comprehensive calibration system to minimize position dependent corrections to 
the reconstructed energy. The spherical central detector will be placed inside a water pool instrumented 
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with PMTs to identify and veto cosmic muons and provide shielding from radioactive backgrounds in 
the underground environment. A muon tracker on top of the detector will further enhance the muon 
identification. 
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Figure 30: Experimental site of JUNO. Hongkong and Daya Bay reactor complex are labeled. Repro¬ 
duced with the permission from Ref. |128j . 
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Figure 31: One conceptual design of JUNO detector. Reproduced with the permission from Ref. |128] . 


RENO-50 is a proposed large (18 kton) liquid scintillator detector in Korea |13Uj . At the present 
time, the schedule for the facility and detector construction aims at years 2015-2020 with the start of 
operation in 2021. The central detector is expected to be equipped with 15k high quantum efficiency 20- 


inch PMTs. Fig. [32] shows the proposed site of RENO-50, roughly 50 km away from the Hanbit nuclear 
power plant. The physics goal of RENO-50 is similar to those of JUNO, including the determination of 
MH and the precision measurements of neutrino mixings. In order to reach the desired energy resolution 
for MH determination, RENO-50 is carrying out R&D to i) enhance the light yield of liquid scintillator 
from ~ 500 to > 1000 photoelectrons/MeV, ii) increase the attenuation length of the liquid scinillator 
from ~ 15 to ~25 m, and iii) enhance the quantum efficiency of PMT from 20% to 35%. Additional 
R&D include the reduction of liquid scintillator radioactivity to 10“^®g/g for U and Th. 
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Figure 32: The RENO-50 detector will be located at underground of Mt. Guemseong in a city of Naju, 
47 km from the Hanbit nuclear power plant. Reproduced with the permission from Ref. |130] . 


Reactor antineutrino spectrum 



Figure 33: The spectral distortion at JUNO (arbitrary scale in the vertical axis) for both normal 
and inverted hierarchies with a perfect energy resolution. Am ‘^2 is assumed to be 2.41x10“^ eV^ and 
—2.53 X 10“^ eV^ for NH and IH, respectively. The reactor spectrum without oscillation and the 
oscillated spectrum at 2 km are shown for comparison. 


The principle of the MH determination in a medium baseline reactor neutrino experiment has been 
illustrated in Sec. Figure ^ shows the distortions of the ideal spectrum at JUNO for both normal 
and inverted MH with the perfect energy resolution. The effective mass-squared splitting Amgg (fast 
oscillation frequency) at low energy (2-4 MeV) is higher than that at high energy (4-8 MeV) for NH, 
and vice versa for IH. As the oscillation at the low energy region (2-4 MeV) is very rapid, a good energy 
resolution is required to extract the oscillation signal. In general, a better than 3%/y^i7(MeV) energy 
resolution has been discussed in the literature as the requirement. In reality, the energy resolution is 
generally parameterized as: 

(37) 


E 


b 


U2' 


Here, the 1/E term represents the contribution from photon statistics; the l/E"^ term represents the 
contribution from PMT dark noise; the constant term represents the contribution from residual non- 
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uniformity after calibration. As the MH signal is around 2.5 MeV, the requirement on the energy 
resolution is naturally converted to better than 1.9% at 2.5 MeV prompt energy. Beside the energy 
resolution, a <^1% absolute energy scale uncertainty [TO tlldlj is also required in order to unambiguously 
extract the MH signal. Such a requirement can be addressed with a dedicated calibration system with 
a small positron accelerator (such as a Pelletron). 
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Figure 34: (Left) The expected prompt energy spectrum of JUNO for six years of data taking with 
a 20 kt detector and 36 GWth reactor power (a total of 100k IBD events). The energy resolution is 
assumed to be about 1.8% at 2.5 MeV. The bottom panel shows the MH signal as the ratio of spectrum 
of inverted hierarchy over the normal hierarchy. (Right) JUNO’s sensitivity evolution with respect to 
calendar years. [73]. The AT is a test statistics consisting likelihoods of normal and inverted MH for 
data X. The green and yellow bands represent the 68% and 95% expectations, respectively, taking into 
account the fluctuations in statistics and variations in systematic. The dashed lines correspond to the 
probability ratios of the normal vs. inverted MH in the Bayesian framework. 


In addition to MH, JUNO will perform precision measurements of the neutrino mixing, which is a 
powerful tool to test the standard 3-flavor neutrino framework (or z/SM). The precision measurement of 
sin^ 2012 will i) lay the foundation for a future sub-1% direct unitarity test of the PMNS matrix [1321 142] . 
ii) constrain the allowed region of the effective neutrino mass associated with the decay rate of the 
neutrinoless double beta decay, and hi) test various phenomenological models of neutrino masses and 
mixing |133j . such as 0i2 = 35°-|-0i3Cos ^12 = 32°-|-0i3 cosJcP) and 0i2 = 45-|-0i3 cosJcP) with 

being the CP phase in PMNS matrix. The precision measurement of Am^g (or Am 32 ) will i) test an 
important sum rule Ami 3 -|-Am^i J-Am |2 ~ 0 and ii) reveal additional information regarding MH, when 
combined with the precision Am^^ measurements from muon neutrino and antineutrino disappearance 
in accelerator experiments. JUNO’s expected precision of Am^i, Am| 2 , and sin^ 202i will be about 
~0. 6% uni. These precision measurements are complementary to the atmospheric and long-baseline 
accelerator neutrino programs. Together, the next-generation reactor neutrino and accelerator neutrino 
experiments will access all seven neutrino mixing parameters governing the oscillation phenomenon 
within the three-neutrino framework. 

To summarize, the major physics goals of the next generation medium baseline reactor experiments 
such as JUNO and RENO-50 are: 
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to be the first experiment to simultaneously observe neutrino oscillations from both the atmo¬ 


spheric and the solar neutrino mass-squared splittings (see Fig. 34), 


to be the hrst experiment to observe more than two oscillation cycles of the atmospheric mass- 


squared splitting (see Fig. 34), 


determination of the neutrino mass hierarchy, i.e. whether is positive or negative, through 


the precision measurement of the spectral distortion. Right panel of Fig. [34] shows the expected 
sensitivity of JUNO vs. the running time. 


precision measurements of sin^26'2i, Am 32 , and Am^i to better than 1%. We should note that 


the precision measurement of Am 32 requires the knowledge of the neutrino mass hierarchy. 


Besides these, the 20 kt underground liquid scintillator detector offers a rich physics program including 
the search for proton decay, geoneutrinos, supernova neutrinos, and many exotic neutrino physics topics. 


4 Summary and Outlook 

In this work we review the theoretical and experimental issues one has to resolve in order to determine 
the neutrino mass hierarchy. The motivations and the theoretical framework are depicted in Sec. [T] 
Then, in Sec. [^ the practical methods that can be used in the MH determination are described. Finally, 
in Sec. [^the various existing and proposed or planned experiments that have the capacity to shed light 
on the MH are briefly reviewed. 

Determining the neutrino mass hierarchy is a challenging problem. Even though one could naively 
expect that the difference between the normal and inverted hierarchies, the two possible mass ordering 
shown in Fig. is going to be signihcant, in fact the difference is rather subtle in the real data. That 
is so because the majority of experimental hndings, at least for the most common case of neutrinos 
propagating in vacuum, can be described in a reasonable approximation as stemming from mixing of 
just two neutrino states; in that case the oscillation pattern depends only on the absolute values of the 
mass-squared differences, not on their sign. Thus, in order to determine the MH, i.e. the sign of Am|^, 
it is necessary to go beyond that approximation and consider the effects associated with mixing of the 
three neutrino states. The corresponding effects are, however, suppressed by the two small parameters 
present in the oscillation formalism, the mixing angle 013 ~ S.b-’ and the mass square difference ratio 
Am|i/Am^2 1/30. The fact that it is possible to approach the problem in this way at all is the result 
of successful experiments that have recently shown that the mixing angle 0 i 3 is not as small as many 
people believed. That Ending opens the way to the MH determination. At the same time, that finding 
also shows that the seemingly simple and attractive possibility of the so-called tribimaximal mixing 
matrix, originally suggested in Ref. ra and discussed in many papers since, is not realized in nature. 
Even though the PMNS matrix, as determined from the oscillation experiments, looks similar to the 
simple form suggested, with |UpP expressed in rational numbers 2/3,1/2, 1/3,1/6 and 0, the reality is 
more complicated. 

The other way to approach the MH is to use the neutrino propagation in matter, with the additional 
interaction that depends on the mass ordering. The mass ordering of neutrinos ui and z /2 (ruj,^ < 
is determined in this way thanks to the results of the solar neutrino experiments. Solar neutrinos 
propagate from the central region of the Sun, over large distance with high electron density, so the 
matter effects are very pronounced. However, the corresponding oscillation length Lmatter to determine 
the mass ordering of neutrinos vi and (see Eq. (|^) is very long O(IO^) km for the typical terrestrial 
electron densities, so the practical application requires very intense neutrino sources and large and 
complicated detectors. 
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Instead of neutrinos from accelerators, one could use the atmospheric neutrinos that come from all 
directions, including the upward going ones that traverse the Earth diameter. In that case the matter 
effect is large, but the MH determination is hindered by the fact that the atmospheric neutrino flux 
contains both and v^. and also both v and v and thus one needs a way to recognize the flavor and 
polarity of the interaction products, or employ the fact that the different components of the atmospheric 
neutrino flux have different intensities and different interaction cross sections. 

Having enumerated the difficulties we can understand why the MH remains unknown. But that 
situation can change. As discussed in Sec. one can approach the MH problem with neutrinos from 
accelerators and from the atmosphere through the matter effect, or from nuclear reactors through 
exploring difference between ^"^31 • Since the effects are subtle, it is necessary to use 

adequate devices and reliable statistical analysis methods. Also, it would be desirable if all, or at 
least several, of the possible methods were used in practice. The results would be complementary and 
increase the confidence that the conclusions are indeed correct. 

While some of the existing experiments, described in Sec. have a chance to resolve the MH problem 
soon, it is likely that eventually the more complex and costly devices described there will be needed. 
Again, complementarity will be desirable and very useful. Realistically, some of the listed experiments 
might begin operation within this decade, with results after few years of operation. So stay tuned. 
Rewards are coming only to those who are patient. 
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